A GENERALIZATION OF STANLEY'S MONSTER RECIPROCITY 

THEOREM 

GUOCE XIN 

Abstract. By studying the reciprocity property of linear Diophantine systems in 
light of Malcev-Neumann series, we present in this paper a new approach to and a 
generalization of Stanley's monster reciprocity theorem. A formula for the "error 
term" is given in the case when the system does not have the reciprocity property. 
We also give a short proof of Stanley's reciprocity theorem for linear homogeneous 
Diophantine systems. 
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1. Introduction 

Let A be an r by n matrix with integer entries, and let b be an r- vector in TIT . Many 
combinatorial problems turn out to be equivalent to finding all nonnegative integral 
(column) vectors a G W 1 satisfying 

(1.1) Aa = b, 

especially in the homogeneous case when b equals 0, of which the solution space is a 
rational cone. Such problems are also known as solving a linear Diophantine system. 
There are two closely related generating functions associated to (jl.lj) : 



£(x;b)= J2 x ? £( x ; b )= E x * 

(ai,...,a„)eN n (ai,...,a„)eP" 



ai - - - x a 



n i 



where the first sum ranges over all a = (a±, . . . ,a n ) such that Aa = b, and the 
second sum ranges over all positive integral a such that Aa = — h. We omit b in the 
homogeneous case. The following well-known reciprocity theorem for homogeneous 
linear diophantine equations was given by Stanley as 0, Theorem 4.1]. 

Theorem 1.1 (Reciprocity Theorem). Let A be an r by n integral matrix of full rank 
r. If there is at least one aeP" such that Aa = 0, then we have as rational functions 

E(x u ... } x n ) = (-l) n - r E(xi 1 ,... 1 x- 1 ). 
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Previous proofs of this theorem used decompositions into simplicial cones or lattice 
cones, or complicated algebraic technique. See 0, p. 214] and p| for further informa- 
tion. We will give a short proof using a signed cone decomposition and induction. 

In the general situation, the best known result (up to now) is the monster reciprocity 
theorem, which was given by Stanley Q in 1974. The theorem will be stated later after 
new notation is introduced. It includes as special cases many combinatorial reciprocity 
theorems, such as the reciprocity theorem for homogeneous linear Diophantine system, 
that for Ehrhart polynomials, and that for P-partitions, etc. We will give a simple 
approach to this theorem. As applications, we give detailed, and short, implication of 
the reciprocal domain theorem j4j, Proposition 8.3]. 

The new approach uses the idea of Malcev- Neumann series 0; 0; 0], which defines 
a total ordering on the group of monomials to clarify the series expansion of rational 
functions. We study the reciprocity property of an object that is more general, but 
less combinatorial, than that was studied in |4j . The new objects we are going to study 
are Elliott-rational functions, while the previous objects are Elliott-rational functions 
with a monomial numerator. By an Elliott-rational function, we mean the one that 
can be written as 



where p is a polynomial and yi and Z{ are monomials. 

In this larger set of objects, it is much easier to build up the reduction steps. Theorem 
13.81 a general result that gives a reciprocity formula for Elliott-rational functions, turns 
out to be easy to prove. We shall use this result to formulate the monster reciprocity 
theorem (Theorem 14.2)1 . 

In Section 2, we introduce the basic idea of Malcev- Neumann series and reformulate 
the reciprocity of linear Diophantine system in terms of constant terms. In Section 3, 
we develop the reciprocity theorem for Elliott-rational functions. We apply our result 
in Section 4 to give the generalized monster reciprocity theorem. In section 5, we 
illustrate the monster reciprocity theorem by examples, and as an application, we give 
a simple derivation of Theorem 11.11 Section 6 includes an inductive (combinatorial) 
proof of Theorem 11.11 



F(Ai,...,A r ,x) 




YlZliVi-Zi) 



2. Reciprocity in Terms of Constant Terms 



Solving a linear Diophantine system (LD-system for short) means finding all vectors 
a G N n that satisfy Aa = b, where A is an r by n matrix with integral entries. More 
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precisely, we want to solve the following system of equations: 

ai,iai + ai j2 a2 + h CL^ n a n = h 

a 2 ,ia 1 + a 2 ,2«2 H h fl2,n«« = ^2 

(2-1) =••• 

Or,l a l + a r,2«2 + h O-r^Otn = K- 

We assume the rank of A\b equals the rank of A, for otherwise, the LD-system has no 
solution even in Q. 

Let Ci be the ith column vector of A. Then the above system is the same as 

Cioci + C 2 a 2 H h C n a n = b. 

Now let E(h) and E(b) be the sets of all such solutions in N n and P n respectively. 
It is interesting to study the following two associated generating functions of (|2.1jl : 

(2.2) £(x; b) = £ x", £(x; b) = ^ x" 

aG-B(b) aeJ(-b) 

where x = (x 1; . . . ,x n ) and if a = («i, . . . ,a n ), then x Q := x" 1 • • -x" n . 

The above equation defines two rational functions in x. If as rational functions 
-E"(x; b) = (— l) n " r £'(x _1 ; b), then we say that the system (|2.1j) has the R-property 
(short for reciprocity property). 

We can compute i?(x; b) by replacing the r linear constraints with r new variables 
Ai, X 2 , . . . ,X r and then take the constant terms. Let A be (Ai, . . . , A r ), and let CTa F 
be the constant term of F in A. We have 



E{x-h) = CTA" 1 ' iai+ '" +ai ' nan " 6l ---A^ iai+ - + ^' ia ' l - b '-x Q 

aeN n 

\— 61 \-b r A-b 

1 j " nr=i(i - *r at • • • a^*,) v rnua - a c ^) ? 

with the working ring C[A, A _1 ][[x]], where A -1 means (A] -1 , . . . , A7 1 ). The above 



conversion can be trait back to MacMahon Similarly we get 

(2-4) E(x;b) = CT n „ 11 ;=^ V c . V 

We define £(A, x; b) and £(A, x; b) to be the crude generating functions of E(x, b) 
and -E(x; b) as 

(2.5) g(A,x;b) = =- 7 - £(A,x;b)- ' 1Ll 



nr =1 a-A^)' v ' ' ^ nr=i(i-A^) 

and observe that as rational functions 

^(A- 1 ,x- 1 ;b) = (-l)^(A,x;b). 
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However, the series expansion of the two sides of the above equation are different. The 
change of variables by A — > A -1 , which corresponds to multiplying each row of (j2.1j) 
by —1, will not make a difference when taking constant terms. Therefore, the system 
has the R-property if and only if as rational functions 

CT£(A,x) = (-l)'CT'5(A,x), 

A A 

where we expand S(A, x) on the LHS at x = 0, while on the RHS at x = oo. 

As we shall see later, the different expansions appearing in the above equation is 
easily explained in the context of Malcev-Neumann series. 

The group of monomials in A and x can be given a total ordering "^ p " that is 
compatible with its group structure; i.e., for any monomials A, B and C, A ^ B 
implies AC ^< p BC. This is equivalent to a total ordering < p on the additive group 
Z n+r . An important such ordering < is the reverse lexicographical ordering on Z n+r . 
Then a Malcev-Neumann series (or MN-series for short) with respect to -< p is a formal 
series on A and x with a well-ordered support: the set of monomials corresponds to the 
nonzero terms. Recall that a well-ordered set is a totally ordered set such that every 
nonempty subset has a minimum. 

For our purpose, p will denote an injective endomorphism of Z ra+r (a nonsingular 
integral matrix), and < p will be the induced total ordering defined by a < p b if and 
only if p(a) < p(b). We denote by C p ((A,x)) the corresponding field of MN-series with 
respect to p. The field of iterated Laurent series C((A, x)), where p is the identity map 
and is omitted, has been studied in (s[ E| . For a more general setting of MN-series, the 
readers are referred to 0; or 3, Chapter 13]. 

The series expansion of MN-series will be explained in more details in the next 
section. Let us review some properties of MN-series Q to see that such fields are 
suitable for dealing with different kinds of series expansions of rational functions. 

For any total ordering < p , C P ((A, x)} is a field. In particular, C((A, x)) is the field of 
iterated Laurent series 

The field C(A, x) of rational functions is naturally embedded into C P ((A, x)) for any 
p. This follows from the field structure of C P ((A, x)) and the fact that every polynomial 
has a finite support. 

Every rational function F(A, x) has a unique expansion in C P ((A, x)). The expansions 
of F for different p are usually different. For instance, the expansion of l/(x — y) in 
K((x, y)) is 



i i 



x — y x 1 — y/x x 



but the expansion in K((y, x)) is 

1 1 



x-y -y 1-x/y -y fc f| 
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Note that we can write K((y, x)) as K p ((x, y)) where p is defined by the matrix 

or by abuse of notation, p(x) = y and p(y) = x. 

Recall also that every subset of a well-ordered set is well-ordered. Thus the following 
operators CTa, PTa, and NT a are well-defined for MN-series. 

CT = b , FT M* = b * Xk > and N A T E hkXk = E hkXK 

k£Z keZ k>0 keZ k<0 

Obviously, for an MN-series F(A), CTa-F(A) = PTa F(A)| A=0 . The constant term 
operators are commutative so that taking the constant term in a set of variables is 
defined by iteration. 

Now it is easy to see that Theorem II .11 is a consequence of the following proposition. 
Proposition 2.1. Suppose that E is nonempty. Then 
(2.6) CT£(x;0) = (-l) rank(A) CT p £(x; 0), 

A A 

where p is the endomorphism defined by p(xi) = x^ 1 and p(Aj) = Aj. 

On the other hand, it is easy to deal with the case of rank(A) < r. So Theorem ll.il 
is equivalent to Proposition 12 .1\ whose proof will be given in section 6. 

The following lemma asserts that elementary row operation will not change the 
solution space of an LD-system. We give it here to show that all the work can be done 
algebraically. 

Lemma 2.2 (0, Corollary 3.18). Suppose y is another set of variables. If $ G 
A^x- 1 ]^)), then for fi = x[ a ■ ■ -x b ™ with det(6 ij ) 1 < iJ < n ^ 0, 

CT$(/i,...,/ n ) = CT$(zi,...,z n ). 

X X 

3. Reciprocity of Elliott-Rational Functions 

It is convenient for our purpose to denote by K the field C(x). The field of rational 
functions C(A,x) can be identified with K(A). Usually we are taking constant terms 
in the ring C[A, A _1 ][[x]]. This ring can be embedded into C p ((A,x)), as long as ^ p 
is compatible with the relation Xi y?~ p Xj for all i and j, where A y?~ p B means that 
A >~ p B k for any positive integer k. 

The case r = 1 is illustrative for our understanding of the series expansion for MN- 
series, and in this particular case, we need not restrict ourselves to Elliott-rational 
functions. Let us consider the following problem. 

Problem: Given a rational function Q(X) (short for Q(A,x)) of A and x, compute 
PT^Q(A,x), where the notation PT^ indicates that Q(X) is treated as an element of 
C P ((A, x)), and we use similar notations for the CT and NT operators. 
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To deal with this problem, we shall understand that Q(X) is not only an element of 
K(X), but also an element of C p ((A,x)}. As an element of K(X), Q(X) can be written 
as p(X)/q(X), where p(X) and q(X) are both in K[X]. As an element of C p ((A,x)), the 
denominator 5(A) plays an important role. 

Recall that C P ((A, x)) is equipped with a total ordering ^ p on its group of monomials 
and that its elements have well-ordered supports. Thus for a nonzero element 77, we 
can define its order ordr] to be minsupp^), and its initial term to be the term with 
the least order. The order of is treated as 00. Let us write 5(A) = J2i=o a i^\ with 
cii G C(x) and 7^ 0. To expand Q(X) into a series in C p ((A,x)), we need to find the 
X-initial term a^A- 7 , i.e., the j such that ord(ajA J ) -< p ord(cijA*) for all i 7^ j. This can 
be achieved because of the different powers in A. Then 

\ \ \ 1 ^ / ^ \ k 

SA) = ^ 1 + Ei#j = ^ g (_1)fe ( g ai/a ^J ■ 

This expansion is justified by the composition law 0, Theorem 2.2]. 
It is now clear that we have the following three situations. 

(1) If j equals 0, then for any polynomial p( A), p(X)/q(X) contains only nonnegative 
powers in A. In this case, we say that l/q(X) is PT P in A. 

(2) If j equals d, then for any polynomial p(X) of degree in A less than d, p(X)/q(X) 
contains only negative powers in A. In this case, we say that l/q(X) is NT P in 
A. 

(3) If j equals neither 0, nor d, then 1/^(A) contains both positive and negative 
powers in A. Thus 1/<?(A) is neither PT P nor NT P in A. 

Lemma 3.1. Let q\ and q2 be polynomials in X. Then for any total ordering < p 

• Both l/gi(A) and l/q2 (A) are PT P in X if and only if 1/(5152) is. 

• Both l/gi(A) and l/q2 (A) are NT P in X if and only if 1/(5152) is- 

• For all the other cases, 1/(5152) is neither PT P in X nor NT P in X. 

Proof. We prove the first case for PT as follows. The other cases are similar. Write 

di d2 di+d-2 

5i = ^a-iX 1 , 52 = ^ hX\ and 5152 = C;A\ 

i=0 i=0 i=0 

Suppose that a^A- 71 and bj 2 X^ 2 are the A- initial term of 51 and 52 respectively. Now if we 
expand the product 5x52 but do not collect terms, then a^&^A- 71 ^ 2 is the unique term 
with the least order. So the order of c^+^A- 71 ^ 2 has to equal the order of a^&^A- 71-1- - 72 . 
This implies that the A- initial term of 5152 is Cj 1+ j 2 X : ' 1+: > 2 . The assertion for PT in the 
lemma hence follows from the fact that j\ + fi — ji — and J2 = 0. (Remember 
that ji,j 2 > 0). □ 

A direct consequence of the above lemma is the following corollary. 
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Corollary 3.2. Ifl/q^X) is PT P in X and l/g 2 (A) NT P A, then qi(X) and q 2 (X) 
cannot have a nontrivial common divisor in K[X], i.e., they are relatively prime. 

Definition 3.3. If g(A) can be factored as gi(A)g 2 (A) such that l/gi(A) is PT P in A 
and l/g 2 (A) is NT P in A, then we say that g(A) is p-factorable, and g(A) = <?i(A)g 2 (A) 
is a p-factorization. Such factorization is unique (if it exists) up to a constant in K. 

Theorem 3.4. Letp(X),q(X) G K[X]. If q(X) is p-factorable, then CT p p(A)/g(A) is 
in K, i.e., is rational. 

Proof. Suppose g(A) = <?i(A)g 2 (A) is such a p-factorization. Since l/gi(A) is PT P in A 
and l/g 2 (A) is NT P in A, gi(A) and g 2 (A) are relatively prime in K[X] by Corollary 13.21 
Thus we have the unique partial fraction expansion in K(X): 

(3.1) MX) + ElW + P^ 

9(A) gi(A) g 2 (A) 

where p, are polynomials in A for i = 0, 1, 2 and degpj(A) < deg g^(A) for i = 1, 2. Since 
when expanded as series in C P ((X, x)), po(X) andpi(A)/gi(A) contains only nonnegative 
powers in A, and p 2 (A)/g 2 (A) contains only negative powers in A, we have 

PT^=p (A) + ^. 
g(A) gx(A) 

Thus CT^ = p (0) +Pi(0)/gi(0) is in C(x). □ 

This theorem generalizes a result of Hadamard [6, Proposition 4.2.5 ], which says 
that the Hadamard product of two rational power series is rational. This statement can 
be easily seen from the following observation: Let / = J2k>of kxk an< ^ 9 = Ylk>o9k xk ■ 
Then the Hadamard product of / and g is 



J2fk9 k x k = CTf(X)g(x/X), 



A 
A:>0 

where we are taking the constant term in C((A, x)) for the RHS of the above equation. 

For any total ordering < p on the monomials of K(X), we let ^ p be the total ordering 
such that mi -< p m 2 if and only if mi >~ p m 2 for all monomials mi and m 2 . 

Then we have a sort of reciprocity invariant, for which we need three more notations. 
We use the notation CT^o-^XA) to indicate that F(X) is treated as an element in 
K((X)) and CT A=00 F(A) to indicate that F(X) is treated as an element in K((A _1 )). 
We define 

(3.2) Xa^(A) = CT F(X) + CT F(X). 

A=0 A=oo 



cS 
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Theorem 3.5. Suppose that p(\),q(X) G K[\], with q(X) being p-factorable. Then 
the following is always true as rational functions in K: 

(3.3) CTP 4^ + CT ^=^- 

a q(\) x g(A) g(A) 

Theorem 13.51 gives an invariant of a rational function when taking the constant 
term in A. This invariant is independent of the choice of the total ordering < p when 
applicable. This fact is the key in our new approach to the monster reciprocity theorem 
in Section 4. 

Proof of Theorem, \S.b\ Write g(A) as gi(A)g2(A)A s , such that l/gi(A) is PT P in A, 
l/g 2 (A) is NT P , and 91(0)^2(0) 7^ 0. Clearly s > and the partial fraction decom- 
position of p(X)/q(X) can be written as 

p(X) _ p-i(A) pi(X) p 2 (X) 

^A)-^ +MA) + ^(A) + ^(A)' 
where degp__i < s, degpi < deggi, degp 2 < degg 2 , and p is a polynomial. 

Now we are going to apply different operators on this partial fraction decomposition. 
Applying CT P top{X)/q{X) gives us Po(0)+pi(0)/gi(0), and applying CT P top(A)/g(A) 
gives us po(0) +p2(0)/?2(0). Therefore 

CTP ka) CT - pP _w = P m + P m 

a g(A) a q(X) F0K } 9l (0) g 2 (0) 
Applying CT A=0 top(A)/g(A) gives us Po(0)+p 1 (0)/q 1 (0)+p 2 (0)/q 2 (0), and applying 
CTa=oo to p(X)/q(X) gives us Po(0). Thus the theorem follows. □ 

Remark 3.6. In the proof of Theorem 13. 5[ we see that p can be replaced with a if a 
switches the PT and NT properties of l/gi(A) and l/q2(X) with respect to p. 

As an element of K[X], q(X) can be factored into the product of irreducible poly- 
nomials. Let q(X) = qi(X) ■ • • qk(X) be such a factorization. By Lemma I3~T1 q(X) is 
p-factorable if and only if every 1/qi is either PT P or NT P . When this is the case, the 
p-factorization can be obtained by collecting similar terms. 

Elliott-rational functions are p-factorable for any p. Such a function F can be written 
as follows: 

(OA) F = 

1 • ) (Aii _ ai ) . . . ( A in _ a n )(A*i - 60 ■ • ■ (X k ™ - b m ) ' 

where p(A) is a polynomial in A, ji and hi are positive integers, m and n are nonnegative 
integers, and a« and b\ are monomials independent of A. For a particular p, we require 
that \ j{X ii - ai) is NT P in A, and l/(A fcl - h) is PT P in A. Note that a x can be 0. "The 
method of Elliott" 0, p. 111-114] shows that CT P F is always Elliott-rational. 

A rational function of A is proper in A if the degree in A of its numerator is less than 
that of its denominator. 
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Corollary 3.7. Let F(X) be of the form (|3.4jl . If F(0) = 0, and F(X) is proper in X, 
then for any p, we have a reciprocity formula 

CT p F{X) = -CT^F(A), 

A A 

where both sides are regarded as elements in K . 

More generally, a rational function F is said to have the R-property with respect to 
pif 

(3.5) CT P F = (-l) d CT^F. 

A A 

for some integer d. Here we restrict our interest to the case when d equals r, the 
number of A's. We have the following reciprocity formula for Elliott-rational functions. 

Theorem 3.8. Let F(A, x) be an Elliott-rational function and let < p be a total or- 
dering on Z n+r that is compatible with its additive group structure. Then 

r-l 

(3.6) CT P F = (-l) r CT p F + y^{-l)' 1 CT P 1 X +1 CT p F, 

A A ^ A r ,...,A i+ 2 Ai,...,Ai 

where CT^. Ai is the identity operator for i = and similar for CT P r A +2 when 
i = r — 1. 

Proof. Since we are always taking constant terms in Aj, we omit the A for convenience. 
We compute the following in two different ways. 

r-l 

(3.7) y"(-l) J CT p CT P F + (— l) 1 CT^ CT P F. 

' ^ r,...,i+l r,...,i+2 A 

i=0 



Using Theorem 13.51 we can rewrite ()3.7|) as 

r-l 

V(-l) i CT^Zi+i CT"F. 

-t— n r„,.,i+2 



i=0 

On the other hand, most of the terms in (|3.7|) cancel with each other. The only 
terms left are given by 

CT P F + (— l) r_1 CT P F. 

r,...,2,l r,...,2,l 

The proposition then follows. □ 

Theorem 13.81 gives the error term of the reciprocity formula. A different error term 
representation was given in p in terms of cohomology. However, the computation of 
this error term saved only a little work for general r. Our formula for the error term is 
true for any fixed order of Ai, . . . , A r , and any fixed order of x%, . . . , x n . This suggests 
that some simplifications might exist and a better formula is possible. We have not 
succeeded in finding a such formula. 
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Since simple equivalent condition for F to have the R-property is unlikely, we search 
for a sufficient condition. Corollary 13.91 and Proposition 13.101 below play important 
roles in our formulating the monster theorem. 

A rational function F is said to have the I-property with respect to p if for i = 
1, 2, . . . , r, we have 

(3.8) l Xt CT p ---CT p F = 0. 

Ai 

Corollary 3.9. If an Elliott-rational function has the I-property, then it has the R- 
property. 

This result is a direct consequence of Theorem 13.81 The special case of this corol- 
lary when the Elliott-rational function has a monomial numerator was shown by a 
complicated computation in 0, Lemma 9.2]. 

In the case r = 1, Theorem 13.81 gives the equivalence between the I-property and 
the R-property. Moreover, we have, as shown below, a nice equivalent condition 0, 
Proposition 10.3] for the R-property that contains no algebraic expression. 

Proposition 3.10 (4]). LetS(x;b) be the crude generating function associated to an 
LD-system consisting of a single equation Aa = a\ct\ + ■ ■ • + a n a n = b: 

^ (X ' 6) = (l- A fl ixi)---(1- A a »x n )' 
Then the following four conditions are equivalent for any p: 

(1) £(x;6) has the R-property. 

(2) £(x;6) has the I-property. 

(3) CT A=0 £(x; b) = and CT A=00 £(x; b) = 0. 

(4) The following two conditions are both satisfied: 

(a) There does not exist a (3 G Z with Af3 = b such that (3 e < if a e > and 
(3 e >0ifa e < 0. 

(b) There does not exist a 7 G Z with A'y = b such that 7 e > if a e > and 
7 e < ifa e < 0. 

The proof of this proposition, which is not given in full detail here, proceeds by 
showing that CT A=0 £(x; b) and CT A=oc £(x; b) have no common terms when expanded 
as Laurent series. The reader is referred to |4j, Proposition 10.3] for details. 

4. The Monster Reciprocity Theorem 

Consider an LD-system Aa = b as in ([2.1| . The crude generating function £(A, x; b) 
is an Elliott-rational function with a monomial numerator. We say such a function has 
the matrix form since we are going to represent it by a matrix. The problem is to find 
a simple sufficient condition for S(A, x; b) to have the R-property. A homology version 
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solution can be found in [5|. The best known result was Stanley's monster reciprocity 
theorem 0, Theorem 10.2], which says that the LD-system has the R-property if certain 
linear combinations of its equations have the R-property. We present here a simple 
approach to this problem. 

The central idea of our approach to this problem, as in 0, is to apply Corollary 
13.91 and Proposition 13.101 If the following checking procedure returns a true, then 
£(A, x; b) has the I-property and hence the R-property with respect to p. Note that 
the converse of this statement is false. 

The checking procedure for S(A, x; b): 

(1) Let Ti = £(A,x;b). If J Al T x (A) ^ then return false. 

(2) Write CT\ 1 P T\ (A) as a sum of matrix forms in an efficient way. For every 
matrix form T 2 , if I\ 2 T 2 ^ 0, then return false. 

(3) Repeat the above step for every matrix form T 2 with respect to A 2 , and then 
for every T 3 with respect to A3, ... , until we have checked if I\ r T r (\ r ) 7^ 0. If 
no false is returned, then return true. 

The basic tool in finding these Tj's is partial fraction decomposition of rational func- 
tions. Using the constant term operators seems neater than using residue operators as 
in 0. 

Our task is to find a simple equivalent condition for the checking procedure to 
return a true. Such a condition will be our monster reciprocity theorem. In order to 
do so, we represent a matrix form T as an augmented matrix. In fact, we can keep 
track everything by adding a row of monomials in the x's on the top and a column of 
monomials in the A's to the left of an LD-system. Therefore, the checking procedure 
will be done by matrix operations. Note that using matrix operations is one important 
aspect of the monster reciprocity theorem. 

We use the following identification: 

2/n+lA- b 



T 



nr=i(i-A^ 



V\ ■■■ Vn 

C\ ■ ■ ■ C n 



Vn+l 
b 



where yi are monomials in x, and Ci are column vectors. It would be clearer if we add 
Aj to the left of the zth row, but this is unnecessary after applying Lemma 12.21 and 
requiring that the ith row (with % > 2) of T s is indexed by A s+i _i. 

The row operations we are going to perform will never involve the top row. The 
column operations, when acting on the first row, are treated as multiplications instead 
of additions for the obvious reason. We alow fractional entries and fractional powers. 
Roots of unity might appear, but will not be a trouble. 

Three special matrix operations will be useful. We define T <— C{i) to be the matrix 
obtained from T by adding — aij/a^j times the ith column to the jth column for all 
j 7^ i. This operation is exactly Gaussian column elimination by taking the (2,i)th 
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entry of T as the pivot. Similarly we define the Gaussian row elimination T <— R(i). 
The third operation T <— D(i) is defined to be the matrix obtained from T by deleting 
the second row and the i-th column. 

Combination of the operations will also be used from left to right. For instance, 
T <— CR(i) := T <— C(i) <— Since row operations commute with column 

operations, we have T <— CR(i) = T <— RC(i). It is easy to verify the following. 



T <- Ci2D(i) = T <- RCD(i) = T <- CD(i). 



For example, if T is given by 



(4.1) 



T 



A x 6 A 2 ' 



(1 - A?xi/A 2 )(l - A 2 x 2 /A 1 )(l - x 3 /A?A 2 ) 





2:2 


x 3 


1 


3 


-1 


-2 


b 


-1 


1 


-1 


c 



then 



T <- C(l) 



Xi a^x-f x 3 xf 



3 

-1 





2 
3 









c + 



Xi x 2 x 3 
3 -1 -2 
I -§ 



1 

c + 



and 



T <— CD(1) =T ^CRD{1) 



rr> rr> 3 ,->-» /-v» 3 

t x J i 2** j \ **j 3 1 



37 1 



c + 



The above three operations are generalized to sequences of integers. For instance, 
T <— R(ii, . . . , i p ) is the matrix obtained from T by applying Gaussian row elimination 
by first taking the (2,ii)th entry of T as the pivot, then taking the (3,i 2 )th entry as 
the pivot, and so on. However, the elimination cannot go backwards. For instance, 
we are not allowed to eliminate the nonzero entries in the second row when taking the 
(3,i 2 )th entry as the pivot. 

More precisely, pick out the i±, . . . , i p ih columns of T±, and rearrange them as follows: 



Ti(h, ...,i P ) ■-- 



x 



»i 



X 



■t'2 



Q2,ii a 2,i 2 
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If all the pivots encountering are nonzero, then when ignoring the top row 



(4.2) Tx(ii, ...,i p )<- R(ii, ...,i p ) = 



will be an upper triangular square matrix followed by a zero matrix, and 



(4.3) . . . , i p ) <- RC(it, ...,i p ) = 



will be a diagonal square matrix followed by a zero matrix. Since the matrix opera- 
tions we have performed do not change the determinants, the a' s is can be inductively 
computed by the formulas a' 1}il = ai,^, and n^=i a 'j,i 3 = ^ e ^( a k,ii)i<k,i<s- 

We denote by M(zi, . . . , Zk) is a generic monomial in z\, . . . , z^ whose exact expres- 
sion is not needed. 

Though we can formulate the monster reciprocity theorem for any p, the result seems 
nicer if we assume that p satisfies the following condition: 

Vy = M(x), yy p \ s ^yM(\ s+1 ,\ s+2 ,...)>-y p Xs, (*)■ 

For example, condition («•) holds for any injective p such that p(xi) is a monomial in x 
and p(Xi) is a monomial in A. We will explain two such p in detail in the next section. 

Definition 4.1. With notation as in (|4.H|) . we define (ii,...,i p ) of distinct entries 
ranging from 1 to n to be a contribution sequence of length p with respect to p if 
y is g s ' ls -< p \ s for all s. The empty sequence is a contribution sequence of length 0. 

The name contribution sequence is in correspondence with the "pole sequence" in 
Q . The condition in this definition will be replaced with simple ones for two special p 
in the next section. 

Theorem 4.2 (Generalized Monster Reciprocity Theorem). Let T be a matrix form 
corresponding to an r by n matrix of full rank, and let ^< p be a total ordering on the 
group of monomials in A and x satisfying condition (*). If for every contribution 
sequence (ii, . . . , i p ) of T with p < r, the second row of T <— RD(ii, . . . , i p ) has the 
R-property, then T has the R-property with respect to p. 



Xi 1 












%%2 






2,2ri 



"p,tp 







Xi 1 


Via ■ 


IJip 


0>l,h 


■ 


■ 







• 





• 


a p,i t 
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Proof of Theorem \4-ty F° r given T and The checking procedure will return a true 
if and only if every T p encountered has the property that 1\ T P = 0, which is the same 
as the condition that the second row of T p has the R-property by Proposition 13.101 

We claim that T p must be similar to the following form for some contribution se- 
quence (ii,..., ip-i): 



(4.4) 



T\ <- CD(i 



) 'p— i 



) = r 1 ^-i2C r D<i 1 ,...,i-_ a ). 



The term similar will be explained later. Assuming T p be given by 1)4.4)1 . we can 
complete the proof of the theorem as follows. We observe that the C operations after 
the R operations do not affect the (p + l)st row (and below) of T\. See (14.2)1 . Therefore 
the second row of T v is the same as the second row of T\ <— RD(ii, . . . , i p -i). 

We prove the claim by induction on p. The claim is trivial for p — 1. Now assume 
the claim is true for p = s and we need to show that the claim is true for p — s + 1. 

By choosing appropriate positive integer N and letting A = \l , (note that T s+1 
will be independt of the choice of N), we can assume that 



t s (\ s ) = t:(\) 



X Vm+l 



n 



X a % 



y-i 

a i 

* 



Vm 
* 



Vm+l 
b 



where a k and b are integers, y k = yk-M(X s+ i, X s+2 , ■ ■ ■ ), and the *'s are column of 
integers that we do not care. Dividing the second row of T' s by iV will give us the the 
second row of T s . Since we have deleted s — 1 columns, m equals n — s + 1. 

We observe that y& = x k 'M(xi 1 , . . . , Xi s _ t ) for k < m, where k' — k equals the number 
of j's such that k' > ij. Therefore yi, ■ ■ ■ ,y m are independent of each other. It is now 
straightforward to check that the partial fraction decomposition of T S '(A) with respect 
to A is: 



T' S (X) = L(X) 



k=l j=l A Sfcj'i/fc 



-CkdVk 1/ak 



~ b y m +i 



where L(X) is a Laurent polynomial in A, and ( k ,j ranges over all a^th roots of unity. 

By Proposition 13.10) I\ a T s = implies that CTa s=0 o T s = and hence CT^ s L(X) = 
0. Together with the fact that for any u independent of A, 



CT"— 

a s A 



if u y p A, 



u 



-u)-\ ifu-^A, 



we have 



\a.k\ 



~-l/o fc N_6~ 

) y-m+i 
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where the sum ranges over all k such that (kjVk -< p A = X^ N , which, by condition 
(*), is equivalent to yf sn< " -< p X s . For such k, we can check that (z'i, . . . , i s -i, k') is 
a contribution sequence. 
Let 

/ ~ — ^-/ a k\ — b ~ 



lit v^(l - Gfc 



where we emphasize Tk as a function of x k }^ ak . By delaying the deletion procedure, 
we can check that 

%(x- 1/ak )=T 1 ^CD(t 1 ,...,i s _ 1 ,k'). 
Then CT> s T s is a sum of T s+1 's, each have the form 

h - — 1 k[.^k,jX k , ) 

for some j with 1 < j < \cik\ and k with . . . , z s _i, k') being a contribution sequence. 

We say that T k is similar to 7^ and it is clear that T k has the I-property if and 
only if Tk has. This completes the proof of the claim. 

□ 

Remark 4.3. If T satisfies the condition in Theorem 14.21 then for p < r, it follows 
from the proof that CT^ Xp T can be expressed as a sum of group terms indexed by 
contribution sequences (zi, . . . ,i p ) of T, with the corresponding group being a sum of 
terms similar to 7\ <— . . . ,i p ). In particular, CT^ can be expressed as at most a 
sum of n(n — 1) • • • (n — r + 1) groups, since there are at most n(n — 1) • • • (n — r + 1) 
contribution sequences. A fast way to compute the sum for each group and an effective 
way to reduce the number of contribution sequences will result in an efficient algorithm 
for computing -E(x, b). 



5. Examples and Applications 

To apply Theorem 14. 21 to a particular LD-system, we need to choose a working field 
C p ((A,x)) to work with. The choice of p is not unique, but we will concentrate on 
two special cases that always work. One is the case that p is the identity; the other is 
equivalent to that in j^. In both cases, we can simplify the condition in finding the 
contribution sequence. 

Case 1: Let C((A,x)) be the working field. The condition y ig s '* a -< \ s in Defini- 
tion 14. II can be replaced with sign(a' s ^^) > 0, where if we write y is = x^ 1 ■ ■ -a^j", then 
/ is the largest such that k\ ^ 0. In practice, we put the sign of k\ at the upper front 
of Vj . 
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Example 5.1. Let (E, (b, c)) the following LD-system: 

3cti — «2 — 2a 3 = 6, 
— oti + a 2 — ctz = c. 

Then the crude generating function T = 7\ of this LD-system is given by (J4.1j) . 
Using Maple, we find that (E, (6, c)) has the R-property for all (6, c) plotted by •, and 
has I-property for all (b, c) plotted by o in the following Figure ^ where we tested all 
— 12 < b, c < 12. Thus the R-property does not implies the I-property. 
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Figure 1. The R-property and I-property for (E, (b, c)). 



Let C((A, x)) be the working field. We want to apply Theorem 14. 21 to find such pairs. 

Since the second row of T\ has only one positive entries, only (1) is a contribution 
sequence of length 1. So after eliminating Ai, we get a sum of three terms similar to 
T 2 given by 



T 2 = T 1 ^ CD(1) 



1 2 

2 _5 

3 3 







.1 i 



_6 

" 3 

6 
" 3 



where we kept the first column to keep track the original column numbers. Now it is 
easy to see that the only contribution sequence of length 2 is (1, 2), though we do not 
need it. 

Therefore, Theorem 14.21 tells us that the LD-system has the R-property if the fol- 
lowing two equations have the R-property: 

3«i — «2 — 2«3 = 6, 

2a 2 — 5«3 = 3c + b. 

Using Maple, we find all such (6, c) as plotted by o in Figure [TJ 



A GENERALIZATION OF STANLEY'S MONSTER RECIPROCITY THEOREM 



17 



Example 5.2. Consider the equivalent LD-system (E'(— c, £>)): 



a>i — a 2 + «3 = — c, 
3«i — «2 — 2«3 = b, 

where we multiplied both sides of the second equation by —1 and switched the two 
equations. 

We need to find (b, c) for S to have the R-property, where 



S 



X\ x 2 x 3 

1 -1 1 

3 -1 -2 



This time we have two contribution sequences of length 1: (1) and (3). Therefore, 
Theorem 14.21 tells us that the LD-system has the R-property if the following three 
equations have the R-property, where the second and third equation are from S <— (1) 
and S <— (3): 

«i — a 2 +«3 = — c, 

2a 2 — 5a 3 = b + 3c, 
5ai — 3a 2 = b — 2c. 

Using Maple, we obtain the same pairs (b, c) as in the previous one, i.e., those plotted by 
o in Figure H All these three equations are needed to apply Theorem l4.2l The following 
coincidence is worth mentioning: if we only consider the second and the third equation, 
we will get all (b, c) plotted by • in Figure [TJ i.e., those (6, c) such that (E', (b,c)) has 
the R-property. Since the first equation comes from the empty contribution sequence, 
we come back to check the previous example, which is obviously not the case. 

Case 2: Let p be the injective homomorphism into C((x, A, £)) by p(xi) = xrf and 
p(\i) = A r _j+i. Then the condition in Definit ion 14 . 1 1 can be replaced with sign(<ia / s f ) > 
0, where d is the total degree of jji s in the x's. Since we only need to keep track of the 
total degree of the x's, the X{ in the top row of T can be replaced with 1. The monster 
reciprocity theorem obtained this way is similar to that of |4|, in which the computation 
used integration along the circles |Aj| = 1 — with 1 >>> e\ 5>> ei 5>> • • • , where » 
means "much greater", and the Xi is taken to satisfy \xi\ = 5 < 1 for some positive real 
number 5. In fact, the condition as in Definition 14.11 was completely written in terms 
of determinants. 

Detailed example for this case, which will not be given here, can be found in 0, 
p. 245]. 

Now let us consider Linear homogeneous Diophantine system (LHD-system for short). 
We shall use Theorem 14 .21 to derive the following theorem, which implies the reciprocal 
domain theorem 0, Proposition 8.3] including Theorem ll.il 



18 



GUOCE XIN 



Theorem 5.3. Suppose T is a matrix form corresponding to an LHD-system of full 
rank. Then for any p satisfying T has the R-property if and only if either CT^ T = 
CT^T = or CT P A T ^ and CT^T ^ 0. 

If we let p be the identity map, then we get Theorem 11.11 If we let p(xi) = Xi for 



1, . . . , p and p{xj) = x i 1 for i = p + 1, . . . n for p wit 



i 1 < p < n, and p{\) = Xi 



for all i, then we will get the reciprocal domain theorem J Proposition 8.3]. 
Proof of Theorem \5.'A If T has the R-property, then 

CT P T = (— l)' r CT P T. 

A A 

The implication thus follows. Now we show the converse is true. 

Obviously we can suppose r > 0, CT^T ^ and CT^T ^ 0. We first show that 
the second row of T has the R-property. Since T corresponds to an LHD-system, we 
can write 

1 



T 



niLi(i-^i i 



where & is a monomial independent of Ai. If some of the a, are positive and some 
of the ^ are negative, then Corollary 13.71 applies and the second row of T has the 
R-property. Otherwise, one of CTa 1= oT and CT\ 1=oa T will be and the other will be 
nonzero. (Note that since the LHD-system has full rank, the case that a* = for all % 
will not happen.) The statement then follows by Proposition 13.101 

Now by Lemma 12.21 if T' is obtained from T by elementary row operations, then 
Qrpp rp t q an( j Qfjp j't o Therefore, the second row of T' has the R-property. This 
means every linear combination of the equations of T has the R-property. Thus the 
theorem follows from Theorem 14.21 □ 

Remark 5.4. The proof of the theorem, together with Remark 14.31 in fact shows the 
following statement: If CT^T ^ and CT^T ^ 0, then CT^ Ap T is proper in all 
Xi for i > p. On the other hand, a simple proof of the statement will lead to a simple 
proof of Theorem ll.il If we restrict ourself in C[A, A _1 ][[x]], the best known proof of 
Theorem 11.11 should be that given by the author in 3, which is included in the next 
section. 

The above remark suggest a way to reduce the number of contribution sequences 
of an LHD-system: Following the notation as in Remark 14.31 since every T p has the 
R-property for X p , we have a choice to choose all those terms with contribution or 
all those terms (with a minus sign) without contribution. The author is managing to 
develop a computer program implementing these techniques. 
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6. Linear Homogeneous Diophantine Systems 

We are concentrating on linear homogeneous Diophantine systems (LHD-systems 
for short), i.e., Aa = 0. Recall that Cj is the ith column vector of A. We omit the 
so that E and E are the sets of all solutions of Aa = in N n and P ra respectively, 
and similar for other notations. Since the proof closely related the linear system and 
its associate generating functions, we restate them as follows. 



(6.1) £(x) = CT£(x), £(x) = (-1)" CTS^ 1 ). 

A A 

We are going to prove Proposition I2.1[ i.e., to show that if E is nonempty, then 

(6.2) CT£(x) = (_i) rank (A) CT P£( X ) ; 

A A 

where we are taking constant term of MN-series and p{x,{) = x^ for all i. 

We shall see that all of the work is done algebraically. First, let us see some facts. 
Exchanging column i and j corresponds to exchanging x, and Xj. Row operations, 
which will not change the solutions of Aa = 0, are equivalent to multiplying A on the 
left by an invertible matrix. This fact can be obtained by applying Lemma | 

Let us see the simple case of r = 1. In this case, £(x) has the form: 



£M 



\ ai Xn 
1 1 



The condition that E is nonempty is equivalent to saying that some of have to be 
positive and some of have to be negative. Thus when written in the normal form of 
a rational function in A, £(x) is proper and its numerator divides A. So Proposition 
12. II follows from Corollary 13.71 

The general case does not seem to work along this line because of two problems. 
One is how to use the conditions that E is nonempty, and the other is how to connect 
to the rank of A. The proof we are going to give uses induction and Elliott's reduction 
identity ^, p. 111-114], which is easy to check and is not given here. 

Clearly if an, . . . , ai jn are all positive or are all negative, then E is empty. So we 
can assume that an > and a\i < 0. Applying Elliott's reduction identity on Ai, we 
get: 

£ ( x ) = i _ ACi+c 2XlX2 (i_ac 1Xi + i _ \c 2X2 ~ X ) I! i _ AC X . 

Now expand £(x) according to the middle term, and denote the resulting three sum- 
mans by £\, £2, and £3 respectively. We have 

(6.3) £ (x) = £1(21, xix 2 , x 3 , • • • ) +£2(^1^2,^2,^3, • • •) - £3(2^2, £3, • • •)• 
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Then these Si are very similar to S. Correspondingly, they are associated to matrices, 
and hence solution spaces that lie in W 1 and F n . More precisely, Si, i = 1,2,3, are 
associated to A 1 = (Ci, C\ + C 2 , C 3 , . . . , C n ), A 2 = (Ci + C 2 , C 2 , C 3 , . . . , C n ), and 
A 3 = (Ci + C 2 , C 3 , . . . , C n ) respectively. Thus Ei,Ei(x) and E^E^x) are defined 
correspondingly. 

Now the matrix A\ is obtained from A by adding the second column to the first; 
the matrix A 2 is obtained from A by adding the first column to the second. They are 
obtained from A through a column operation. So the rank of A\ and A 2 are both equal 
to that of A. The rank of A 3 might not equal the rank of A. 

Applying CTa and (— l) n CT^ to ()6.3|) respectively, we get our key induction equa- 
tions. 

(6.4) E(x) = E 1 (x 1 ,x 1 x 2 ,x 3 , ...) + E 2 (x 1 x 2 ,x 2 ,x 3 , . ..) - E 3 (x 1 x 2 ,x 3 , . . .), 
E(x) = E x (xi,x x x 2 ,x 3 , ...) + E 2 (x 1 x 2 ,x 2 ,x 3 , . . . ) 

(6.5) + (-i) rank ^- rank ^£ 3 (*iX2, x 3 , . . . ). 

Looking more closely at these Ei, we can see that up to isomorphism, E\, E 2 , 
and E 3 are obtained from E by intersecting the half spaces ct\ > a 2 , a\ < a 2 , and 
the hyperplane a± = a 2 respectively. For instance, (a%, a 2 , ■ ■ ■ , ) belongs to E with 
«i > ot 2 if and only if (a± — a 2 , a 2 , . . . ) belongs to E\. Thus Elliott's reduction identity 
in fact corresponds to a signed decomposition of E. Equation ()6.4|) and ()6.5j) could be 
explained directly from geometry. 

We need two more lemmas to give our proof of Proposition 12.11 We shall see that 
the condition on E plays an important role. 

If E is nonempty, then dimE" = dimi? = n — rank(A). Clearly, the dimension of 
the solution space of Aa = is n — rank(A). Let 7 £ E, and let Ti, . . . , T n _ ran k(A) 
be a Z-basis of the solution space in Z n with Ti = 7. Then for sufficiently large m, 
7717 + Ti, . . . , m7 + T n _ rank ( J 4) will be a linearly independent set in E. 

Lemma 6.1. Suppose that E is nonempty, and that Ei is defined as above for i = 
1,2,3. Then any two of the Ei being nonempty implies that they are all nonempty. 

Proof. Suppose that E\ and E 2 are nonempty. Then we have elements f3 and 7 in E 
such that (3 = (fti, /3 2 , . . . ) with (3i > f3 2 and 7 = (71,72,...) with 71 < j 2 . Then 
(72 — 7i)/3 + (Pi — ^2)7 is in E with the first two entries being equal. This means E 3 
is nonempty. 

Suppose that E\ and E 3 are nonempty. Then we have elements P and 5 in E such 
that P = (Pi, p 2 , . . . ) with Pi > p 2 and 5 = (5i,5 2 ,...) with 8\ = 5 2 . Then for 
sufficiently large m, m5 — P is in E with the first entry being smaller than the second. 
This means E 2 is nonempty. 
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The case that Ei and E 3 are nonempty is similar to the previous case. □ 

Lemma 6.2. If all of the Ei are nonempty, then rank(A3) = rank(v4). 

Proof. By hypothesis, it is clear that E is not contained in the hyperplane oi\ = ol%. 
Thus the intersection of E with the hyperplane has dimension dimE — 1. So dimE 3 
is also dimE — 1 and the rank of A 3 equals n — 1 — dimE^ = rank(A). □ 

Proof of Proposition Wl[ The base case, when A is the zero matrix, is trivial. 

By exchanging rows, we can assume that not all of the entries in the first row of A 
are zero. Moreover, since the entries can not be all positive or negative, we can assume 
the first entry is positive and the second is negative by exchanging columns. 

We use induction on Si(A), which is defined to be the sum of the absolute values 
of all the entries in the first row. Now the above argument applies, and it is easy to 
see that Si(Ai) < Si(A) for i = 1,2, 3. Applying Lemma lQ| we can reduce the seven 
cases of Ei being nonempty or not into the following four cases: 

Case 1: only E\ is nonempty. Let f3 in E be such that (3\ > $2- We claim that 
all a with Aa = satisfy the condition ot\ > a 2 , so that E 2 {x\X2, £2, £3, . . . ) equals 
E s (xiX2, X3, . . .), and hence by induction we have 

E(x) = E 1 (x 1 ,x 1 x 2 ,x 3 ,. ..) 

= (-i) rank(n - Ai) E 1 (xr 1 ,xr 1 x^ 1 ,x 3 : 1 ,...) = (-i) n - rank ^)s(x- 1 ). 

If the claim does not hold, then a>i < a 2 . But for sufficiently large m, m(3 — a will 
produce an element in E 2 or E 3 , a contradiction. 

Case 2: only E2 is nonempty. This is similar to case 1. 

Case 3: only E3 is nonempty. This means that E is contained in the hyperplane 
ai = a 2 . Thus 

E 1 (x 1 ,x 1 x 2 , x 3 ,...) = E 2 (x 1 x 2 , x 2 , x 3 ,...) = ^3(^1^2, x 3 , . . . ), 
and we have 

rank(A 3 ) = n — 1 — dim(E 3 ) = n — dim(E) — 1 = rank(A) — 1. 

So 

E(x) = E 3 (x 1 x 2 ,x 3 , . . .) 

= {-ly^-^^E^x^x^iXz 1 , ...) = (-l) n - rank ( A )£(x- 1 ). 
Case 4: all of Ei are nonempty. By induction, we see that 

Ei(x) = (-l) n - raDk ^)s i (x- 1 ) 

for i = l,2, and that 

E 3 (x 2 , x 3 ,...) = {-\) n - l -^ A *E{x2\ x 3 \ ...). 
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From Lemma \ti.2\ rankf^) = rank(A). Thus together with our key induction equa- 
tions JE3J and (f6~5j) . we get 

£(x) =Ei(xi, xix 2 , x 3 , . . . ) + E 2 (x 1 x 2 , x 2 , x 3 ,...)- E 3 (x 1 x 2 , x 3 ,...) 

=(-i)™ k(A) (^i(xr 1 ) xr 1 x 2 - 1 ,x 3 - 1 , . . .) 

+E 2 (x^ 1 X2 1 , X2 \ xj 1 , . . . ) + i^i" 1 ^ \ ^3 \ • • • )) 

= (-l) n - rankA E(x). 

□ 

Acknowledgement: I am very grateful to Richard Stanley for introducing me to 
his inspiring work. 
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